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ABSTRACT 

Context. Theory suggests that about 10% of 5vyiy/-detected gamma-ray bursts (GRBs) will originate at redshifts, z, greater than 5 yet 
a number of high redshift candidates may be left unconfirmed due to the lack of measured redshifts. 

Aims. Here we introduce our code, GRBz, a method of simultaneous multi-parameter fitting of GRB afterglow optical and near 
infrared, spectral energy distributions. It allows for early determinations of the photometric redshift, spectral index and host extinction 
to be made. 

Methods. We assume that GRB afterglow spectra are well represented by a power-law decay and model the effects of absorption due 
to the Lyman forest and host extinction. We use a genetic algorithm-based routine to simultaneously fit the parameters of interest, and 
a Monte Carlo error analysis. 

Results. We use GRBs of previously determined spectroscopic redshifts to prove our method, while also introducing new near infrared 
data of GRB 990510 which further constrains the value of the host extinction. 

Conclusions. Our method is effective in estimating the photometric redshift of GRBs, relatively unbiased by assumptions of the 
afterglow spectral index or the host galaxy extinction. Monte Carlo error analysis is required as the method of error estimate based on 
the optimum population of the genetic algorithm underestimates errors significantly. 
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1. Introduction 

Theory suggests and observations agree that approximately 10% 
of 5w//f-detected gamma-ray bursts (GRBs) will originate at red- 
shifts, z > 5 (Bromm & Loeb 2006; Jakobsson et al. 2006b) and 
about 3% at redshifts, z > 6 (Daigne, Rossi, & Mochkovitch 
2006). Furthermore, due to the favourable /iT-coiTection at a fixed 
observer time and cosmological time dilation, the observed op- 
tical flux is not expected to fade significantly with increasing 
redshift. Hence GRBs should be detectable out to z > 10 (Ciardi 
& Loeb 2000; Lamb & Reichart 2000), making them the most 
distant observable objects in the Universe. Yet a number of high 
redshift candidates may be left unconfirmed due to the lack 
of observed counterparts or lack of observations (Ruiz-Velasco 
et al. 2007). 

The most usual way of determining the redshift of a GRB 
is to spectroscopically measure the redshift of the afterglow at 
early times if it is bright enough, or to measure that of its host 
galaxy, either photometrically (e.g. Bolzonella, Miralles, & Pello 
2000) or spectroscopically, once the optical afterglow of the 
burst has dimmed sufficiently. Neither of these may be possi- 
ble given the dimness of the afterglows (Fynbo et al. 2007) and 
the host galaxies. In addition, programs to obtain such measure- 
ments are frequently only triggered if there is evidence to sug- 
gest a high redshift, or a burst of significant interest. There are a 
number of redshift indicators, or pseudo redshift methods, which 
rely on correlations between observable properties and redshift 
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(e.g., Guidorzi 2005; Amati 2006; Pelangeon & Atteia 2006) that 
offer approximate redshifts to varying degrees of success. 

An alternative method of redshift determination, useful when 
the source is not bright enough for spectroscopic observations, 
is to photometrically measure the redshift of the afterglow itself, 
given enough simultaneous optical/near infrared data points. 
Attempts at photometric redshift determination usually assume 
a spectral index and a value for host extinction, while fitting for 
redshift, or fit for a limited number of discrete values of spectral 
index and host extinction (e.g., Andersen et al. 2000; Jakobsson 
et al. 2006b). Independently, many attempts have been made 
to fit host extinction and spectral index to similar data, when 
the redshift is previously known (e.g., Galama & Wijers 2001; 
Stratta et al. 2004; Kann, Klose, & Zeh 2006; Schady et al. 2007; 
Starling et al. 2007). This is important as the properties of the 
circumburst medium and host galaxy can be used to constrain 
progenitor types and burst models themselves. 

In this paper we present our code, GRBz, which, unlike other 
methods of photometric redshift determination, fits all three pa- 
rameters simultaneously and allows for determinations of the 
photometric redshift to be made, relatively unbiased by assump- 
tions of the spectral index or the host extinction. In section 2 we 
describe the method of modelling and fitting that we apply to 
the data described in section 3, including our previously unpub- 
lished near IR data of GRB 990510. In section 4 we present our 
results and in section 5 we discuss the results as they apply to 
the fitting mechanism. We summarise our findings in section 6. 
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2. Method 

2.1. Model 

Our method is based on calculating the flux at a frequency, F(v), 
for a given set of fit parameters: redshift, z, spectral index of af- 
terglow, p and host extinction, Eb-v- This is done by assuming 
that the emitted flux from the afterglow is well represented by a 
power-law decay {F(v) oc v^^), and adjusting this for extinction 
in the host galaxy, line blanketing associated with the Lyman for- 
est in the intervening medium and extinction in the Milky Way. 
This flux is then integrated, numerically, over the transmission 
curve of an individual filter, T'fiiter(^): ^min ^ v < Vmax) 

and nor- 
malised to give the flux as measured by that filter. 



Table 1. Optical SED references for the GRBs in our sample. 



^'filter 



1 F(v)T fiitei(v)dv 



Av 



where the efl'ective width of the filter is 

TmteT(v)dv. 

This is repeated for each of the filters in the Spectral Energy 
Distribution (SED) so that the calculated fluxes may be com- 
pared to those measured, by a Chi Squared (x^) test, and best fit 
parameters estimated. 

The extinction in the host galaxy and the Milky Way are cal- 
culated using the models of Pei (1992) or Calzetti et al. (2000). 
We use the Pei model for the MW, SMC or LMC, or the Calzetti 
model for a star forming region, with the corresponding values of 
Rv = Av/Eb-v (CardeUi, Clayton, & Mathis 1989) from those 
authors, as an approximation of a GRB host galaxy. Note that 
here we use Eb-v to parametrise the extinction - not Ay as is 
often used - as it follows naturally from Pel's treatment and sim- 
plifies the calculation of extinction at a given frequency. Ay. 

The Lyman absorption, by neutral hydrogen in the inter- 
galactic medium, is calculated from the model presented by 
Madau (1995). Though Madau only gives the first 4 coefficients, 
Aj (Lyman a,fi, y, S) of the model, the remaining 13 (e...p) may 
be extrapolated from those in Madau et al. ( 1 996) due to their lin- 
ear relationship with the wavelengths from the Balmer formula 
for the Lyman series: 

\-H 

where n = 2, 3, ... , 18 and R is the Rydberg constant. This does 
not take into account the effect of a damped Lyman alpha ab- 
sorber, which may or may not be present in the host galaxy. If 
present this would cause the flux to be under corrected and hence 
cause an overestimate of the redshift. 

As the transmission curves of the filters are read into the 
fitting program directly and the necessary effective frequencies 
and widths calculated, any filter can be used. We have initially 
used flie Johnson UBVRIK (Johnson 1965), 2MASS JHKs 
(Cohen, Wheaton, & Megeath 2003) and Swift UVOT filters 
(Roming et al. 2005), covering the approximate frequency range 
10'"^- 10'^ Hz. 

Independently of the fitting mechanism, it is useful to be 
able to approximate the relative efi'ect that Lyman absorption, 
7'Lyman(>', z), wiU have on observed flux; our code also allows us 
to calculate an effective transmission of a filter at a given red- 
shift, z, and for a given spectral index, j8: 



GRB 


Filters 


Reference 


990510 


Ks,H,J,I,R,V,B 


1 


000131 


K,H,1,R,V 


2 


050319 


V, B, U, UVWl, UVM2, UVW2 


3 


050814 


K,J,I,R,V 


4 


050904 


K,H,J,I,R,V 


5 



1 Stanek et al. (1999); Holland et al. (2000); this article, 2 Andersen 
et al. (2000), 3 Mason et al. (2006), 4 Jakobsson et al. (2006b), 5 
Tagliaferri et al. (2005) 



-^effectiveLyman(z> Py ^filter) 



r°'" rfiiter(v)rLyman(v,z)v '^dv 

^Vmin 

C"" Tatei(v)v~'^dv 



This gives a correction factor by which to divide the observed 
flux so as to estimate the flux as if unaff'ected by Lyman absorp- 
tion, as successfully implemented in Starling et al. (2007). 



2.2. Fitting 

The fitting method implemented within our C-program is 
the genetic algorithm-based optimisation subroutine PIKAIA 
(Charbonneau 1995), which we use to minimise the of the 
data points. Due to the complexity, and many local minima in our 
solution space, this proved to be a much more robust, stable and 
reproducible method than those based on simulated annealing 
or downhill simplex methods (section 10.9 of Press et al. 1992 
and references therein and Nelder & Mead 1965, respectively). 
This method allows us to fix, or constrain, the possible values 
of solutions - for example, if a spectroscopic redshift has been 
determined or the spectral index has been constrained from a 
temporal index - while not being overly dependent on the start- 
ing values. (For an example of the previous use of the genetic 
algorithm in astronomy see Mokiem et al. 2005). 

GRB afterglows can be faint and fade rapidly, so often only 
upper limits on the afterglow brightness are available and we 
require a standard procedure for dealing with this. In these cases 
we take the flux to be zero and use the limit as an estimate of 
error. We have compared this method to a method where thex^ 
is defined as zero when the model flux is within the limit and 
very large elsewhere, and find that the methods give consistent 
results. 

Uncertainties of the fit parameters are estimated via a Monte 
Carlo analysis, whereby the data are randomly perturbed within 
the Gaussian distribution of their errors, and refit multiple times. 
Due to the slowness of these calculations caused by the need to 
numerically integrate over the filters multiple times we choose a 
modest number (~ 10^) of trials. The average of these perturbed 
fits should agree with the initial unperturbed fit and the distribu- 
tion of the parameters should give a accurate approximation of 
the uncertainties. From the frequency distribution of any one pa- 
rameter, we see that the distribution is clearly not Gaussian so the 
usual method of quoting sigma errors is not valid, though we do 
quote the uncertainties as nominal, symmetric 68% confidence 
intervals. A population analysis method may also be used to 
approximate errors directly from the fitting algorithm (Mokiem 
et al. 2005) but the Monte Carlo method is more robust. See also 
Figure 5 and section 5.2. 



p. A. Cuiran et al.: Robust photometric redshift determinations of GRB afterglows at z ^ 2 



3 



Table 2. Near IR Observations of GRB 990510. The magnitudes are 
uncorrected for the Galactic extinction, Eb-v = 0.203. 





Seeing 


Band 




Magnitude 


days 


It 




X sec 




0.601 


0.68 


J 


10x60 


17.51 ± 0.07 


1.000 


1.00 


J 


10x60 


18.16 ± 0.08 


3.837 


1.44 


J 


26x60 


20.66 ± 0.24 


0.610 


0.62 


H 


10x60 


16.96 ± 0.08 


1.014 


1.26 


H 


10x60 


17.52 ± 0.09 


3.862 


1.32 


H 


20x60 


19.92 + 0.28 


11.815 


0.65 


H 


60x60 


22.10 ±0.34 


0.618 


0.63 


Ks 


10x60 


16.25+0.17 


1.024 


1.06 


Ks 


10x60 


16.85 ±0.18 


3.881 


1.34 


Ks 


20x60 


19.25 ± 0.25 



3. Data 

3.1. Sample selection 

The purpose of this paper is to introduce our method of photo- 
metric redshift determination, not to provide a statistical analysis 
of a large number of bursts, so we have "cherry-picked" a small 
sample of bursts (Table 1). GRBs were selected from the litera- 
ture on the basis of availability of an SED and, knowing that the 
Lyman forest would only affect optical bands at z > 2.5, redshift. 
Bursts with previously determined properties, especially a spec- 
troscopic redshift, that could be compared to our fit parameters 
were preferred though in the case of GRB 050814 only a photo- 
metric redshift was available. The previously published optical 
data for GRB 9905 10 were augmented by our near infrared (nIR) 
observations of the source. 

3.2. Near IR observations of GRB 990510 

From May 10 to May 22, 1999, at 4 epochs, a total of 186 
60-second exposures, or frames, in J, H and Ks filters, were 
obtained of the optical afterglow of GRB 9905 10. The data 
were obtained with the Son of ISAAC (SofI) infrared spectro- 
graph and imaging camera on the 3.58m ESO-New Technology 
Telescope (NTT). The NTT-Sofl data were reduced using the 
IRAF package wherein flatfielding, sky subtraction and frame 
addition were carried out. Relative PSF photometry was carried 
out on the final images using the DAOPHOT package (Stetson 
1987) within IRAF. The PSF model was created using 16 stars in 
the field, one of which was a 2MASS object suitable for calibra- 
tion (2MASS designation: 13380057-80291 19, RA Dec (J2000): 
13:38:00.58 -80:29:11.9). The resultant magnitudes and la er- 
rors of the nIR counterpart are shown in Table 2. 

There is clearly a break to a steeper slope between our sec- 
ond and third epochs of observation (Figure 1). This is consis- 
tent with the measurement of a jet break at 1.31 days after the 
trigger time (Zeh, Klose, & Kann 2006). Using the simultane- 
ous fit method outlined in Curran et al. (2007) and this break 
time we find that the temporal indices are ai = 1 .06 + 0.20 and 
a2 = 1.82 ± 0.18 (Icr uncertainties). These values differ from 
those found by Zeh, Klose, & Kann (2006), most Ukely because 
of the lack of temporal sampling of the nIR light curve. 

The nlR magnitudes of GRB 990510 were interpolated to a 
common time (f = 0.61 days) and then converted to flux values 
using the caUbration of Cohen, Wheaton, & Megeath (2003). 
VRI data at the same epoch were calculated from the light curve 
fit of Holland et al. (2000), while B data was calculated from the 
light curve of Stanek et al. (1999). All data were then corrected 




0.5 I 2 5 10 



Time (days) 

Fig. 1. Near IR light curve of GRB 990510 with the simultaneous fit 
shown. 

for Galactic extinction of Eb-v = 0.203 (Schlegel, Finkbeiner, 
& Davis 1998) at their various frequencies (Pei 1992). 

4. Results 

For all the fits presented here we have assumed that the host 
galaxy extinction is well represented by the Pei (1992) SMC 
model. This is supported by studies of GRB host extinction 
(Galama & Wijers 2001; Stratta et al. 2004; Kann, Klose, & Zeh 
2006; Schady et al. 2007; Starling et al. 2007), though it is not 
necessarily the model which best fits each of these individual 
bursts (e.g., GRB 050904; Stratta et al. 2007). As explained in 
section 2.2, best fit parameters are given as the average values of 
those returned from a Monte Carlo analysis (Table 3), uncertain- 
ties are quoted as symmetric 68% confidence intervals. Figure 2 
shows the SEDs for each of the bursts in our sample. The flux is 
calculated from the best fit parameters and clearly shows the cut 
ofl" due to the line blanketing associated with the Lyman forest. 

4.1. GRB 990510 

GRB 9905 10 has a known spectroscopic redshift, Zspec = 1.619± 
0.002 (Vreeswijk et al. 1999; Vreeswijk et al. 2001) which is be- 
low the operational limits of our program, hence we constrain the 
redshift in our fit to be within the errors of the known value while 
leaving the other parameters free. We find that yS = 0.35 + 0.17 
and Eb-v < 0.05, corresponding to a rest frame extinction, 
Av < 0. 14. This spectral index is lower than, but consistent with, 
the value found by Starling et al. (2007) of 0.53+q who utilised 
X-ray measurements as well as our nIR data. Our estimate of 
host extinction is consistent with that of Starling et al. and Stratta 
et al. (2004) who also utiUsed X-ray measurements; all of which 
indicate negligible extinction, though Starling et al.'s approxi- 
mation is much more tightly constrained than ours. 

4.2. GRB 0001 31 

While GRB 000131 has a spectroscopic redshift of Zspec = 
4.500 + 0.015 (Andersen et al. 2000), we have left it, as well 
as the other parameters in our fit, free. From Andersen et al.'s 
published SED, we have found that z = 4.2 ± 0.4, /3 = 0.7 + 0.4 
and Eb-v = 0.10 ± 0.06 (Ay = 0.29 ± 0.18). This extinction 
is consistent with that found by Andersen et al. who assumed a 
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Fig. 2. Fits to the SEDs of our sample GRBs (for filters used in each 
SED see Table 1). The errors on the frequency axis represent the width 
at Transmission = 0.2. 



spectral slope, j8 = 0.70 based on the constraints imposed by the 
temporal decay index. 



4.3. GRB 05031 9 

The UVOT telescope on Swift observed this burst in the UV and 
optical at early times (ro + 240-290 s). Using the pubhshed SED 




10'" 



2x10"' 



5x10'" 
Frequency (Hz) 



lO'S 



Fig. 2. Continued 



(Mason et al. 2006) we find that z = 3.0 ± 0.2, yS = 0.6 ± 0.3 and 
Eb-v < 0.07 (Ay < 0.20). Jakobsson et al. (2006a) determine a 
spectroscopic redshift of Zspec = 3.24 consistent with our result, 
while Mason et al. (2006) infer a spectral index of p - 0.8 + 0.1 
from the X-ray light curve. Though our spectral index is consis- 
tent with this, we cannot infer much information about the spec- 
tral index or extinction since there is only one data point above 
the Lyman break. 



4.4. GRB 05081 4 

Jakobsson et al. (2006b) find a photometric redshift, z - 5.3 ± 
0.3, and fixing yS = 1.0 find an extinction of Ay = 0.9. Using 
their data points we find a redshift of z = 5.77 ± 0.12. Since 
there are only 2 data points above the break, the spectral index 
and extinction may only be loosely constrained as /? = 0.9 ± 0.5 
and Eb-v = 0.08 ± 0.05 (Ay = 0.23). This spectral index is 
consistent with that of the X-ray index, fix = 0.8 ± 0.2 (Morris 
et al. 2005). Jakobsson et al. point out that their extinction is 
marginally higher than inferred from other bursts with bright 
optical counterparts (Kann, Klose, & Zeh 2006). While a low 
extinction is not necessarily the case, our result is in Une with 
those that Jakobsson et al. compare theirs to. These authors also 
point out that their extinction and spectral index would overesti- 
mate, shghtly, the X-ray flux at the time. The alternative spectral 
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GRB 


^literature 


z 


Eb-v 


P 


990510 


1.619±0.002'' 




<0.05 


0.35 ±0.17 


000131 


4.500±0.015'' 


4.2 + 0.4 


0.10 + 0.06 


0.7 ± 0.4 


050319 


3.240±0.001'^ 


3.0 + 0.2 


<0.07 


0.6 ± 0.3 


050814 


5.3±0.3'' 


5.77 ± 0.12 


0.08 ± 0.05 


0.9 ± 0.5 


050904 


6.30±0.002^ 


6.61 ± 0.14 


<0.03 


0.05 ± 0.02 



Vreeswijk et al. (1999), " Andersen et al. (2000), " Jakobsson et al. (2006a), ^ Jakobsson et al. (2006b), ^ Kawai et al. (2006). 
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Fig. 3. Spectral index, j3, versus host extinction, Eb-v, for the Monte 
Carlo analysis of GRB 050814. 

index they suggest to compensate for this (fi = 1.1) is within the 
uncertainties of our estimate. 

4.5. GRB 050904 

The highest redshift GRB identified to date has been the spec- 
troscopically measured, z = 6.295 + 0.002 (Kawai et al. 2006), 
GRB 050904. Using the data pubhshed by Tagliaferri et al. 
(2005), we find a redshift of z = 6.61 +0.14 which is consis- 
tent at the 2cr level with the spectroscopic value. Since the nIR 
flux points, unaff"ected by the break, are all consistent with a flat 
spectrum we can infer very little from the obtained values of 
jS = 0.05 + 0.02 and Eg-v < 0.03. 

5. Discussion 

5.1. Limitations 

We are able to draw a number of general conclusions from the 
data sets that we have analysed here, together with synthesised 
data sets that we initially tested the fitting method on, regard- 
ing the operational limits of the program. Obviously, whether or 
not any one of the three parameters with physical meaning (z, 
P and Eb-v) can be fit with accuracy is dependent on the data 
available. We may only fit the redshifts of bursts with z ^ 2.5 
(lower if UVOT data are available) as this is where the Lyman 
forest starts to affect the high frequency filters (UB). Redshifts 
lower than this act degenerately, so a fit that has a wide spread of 
redshifts under, or around 2.5 must be assumed to be under this 
limit and can only be fit, in the other parameters, if the redshift 
is previously known (e.g. GRB 990510). 



The accurate determination of /3, Eg-v or host extinction 
model is dependent on the number of data points above the 
Lyman break. If there are few data points unaffected by Lyman 
absorption we cannot expect to draw any firm conclusions re- 
garding /3 or Eb-v as the information is lost blue-ward of the 
break. Likewise there must be data points, or limits, blue-ward 
of the break so as to constrain redshift. 

Plotting the various fit parameters for individual bursts 
against each other we see that there may be a dependency be- 
tween host extinction and spectral index (as an example see 
GRB 050814, Figure 3): increasing the spectral index of the 
source, reduces the need for host extinction at a given redshift. 
There may also be a dependency between the redshift and host 
extinction, as a high redshift shifts highly extincted light into the 
observed range, though this efl'ect is not obvious in our results. 

These dependencies should be kept in mind when using the 
best fit parameters but because in our Monte Carlo error analysis 
no one parameter is fixed, as it would be in a normal Ax^ anal- 
ysis, these dependencies do not affect the quoted uncertainties. 
From Figure 3 one can also see the hard limits set on the value 
of the spectral index in the fitting routine for this particular SED: 
0.1 < P < 1.7. Limits are set, out of computational necessity, to 
limit the parameter space being searched by the fitting routine. 
These Umits must be chosen carefuUy, after an initial trial fit, so 
that they include a realistic range of parameters; wide enough to 
include the best fit parameters of each Monte Carlo trial but no 
so wide that the best fit parameters will be difficult to localise in 
the parameter space. 

5.2. Errors 

In genetic algorithm-based optimisation, the population refers to 
every set of parameters that were fit in the search for the mini- 
mum though in the genetic algorithm it is 'fitness' (inversely 
proportional to x^) that is maximised. By plotting the distribu- 
tion of the fitnesses of the entire population, one can see that 
there are many fits of poor fitness, and a number of high fitness 
that returned the best fit parameters. By filtering on that subset 
of the population that were most fit (the optimum population), 
one finds a distribution of the interesting parameters, the width 
of which is taken as an estimate of the uncertainty of the fitting 
mechanism. 

To illustrate, we show these distributions for one of the above 
bursts, GRB 050814 (Figure4). We, arbitrarily though conserva- 
tively, take the optimum population to be those with a fitness 
greater than 0.8. This optimum population consists of two sub- 
populations with fitnesses centred on ~ 0.9 and ~ 1 .0, and these 
two populations can also be seen in the distribution plots of red- 
shift and host extinction. From these plots, we can estimate the 
best fit parameters of GRB 050814 to be: z ~ 5.744 + 0.003, 
Eb-v ~ 0.067 + 0.001, ~ 1.1434 + 0.0002. The optimum 
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Redshift. z 
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Host extinction, Eg_y 



1.1432 1.1434 1.1436 1.1438 

Spectrol index, /? 



Fig. 4. Distribution plots for the initial genetic algorithm fit of GRB 050814. Redshift, host extinction and spectral index distributions are plotted 
for members of the optimum population with a fitness greater than 0.8. 



population width method of uncertainty estimation as described 
by Mokiem et al. (2005) is significantly quicker than a Monte 
Carlo analysis as it requires only one trial of the fitting mech- 
anism, however, if compared to the parameters for this burst in 
Table 3, the errors are also significantly underestimated. They 
are more likely an estimate of the accuracy of the fitting mech- 
anism in finding the nominal minimum for a given set of data 
points, than an estimate of the uncertainties due to the errors on 
the data points. 

The difference in error estimates between the optimum pop- 
ulation width method and the Monte Carlo analysis is clearly 
visible in Figure 5 where we plot the distributions for all trials of 
our Monte Carlo error analysis of this burst, overlaid on the re- 
binned distributions of the optimum population of the first trial 
(Figure 4). Though many of these trials return unacceptable 
(plotted as opposed to the fitness of the population width), they 
are the best fits for the given data and the average parameters 
and errors do not show any deviation when the unacceptable tri- 
als are removed. The distributions of the interesting parameters 
are clearly not Gaussian and hence the usual method of quoting 
sigma errors is not valid. We do quote the uncertainties as sym- 
metric 68% confidence intervals but it is more correct to judge 
the errors from these distribution plots. One can see that the host 
extinction and spectral index of this particular burst have almost 
flat distributions so it would be more proper to give each as a 



range as opposed to central values with errors: Eb-v ^ 0.17, 
0.1 < P < 1.7 (which corresponds to the hard limits imposed 
on jS). Likewise the redshift would be better described with non- 
synmietric, though still non-Gaussian, errors as ~ 

We find similar situations, regarding the underestimate of 
errors by the optimum population width method and the non- 
Gaussian Monte Carlo errors, for the other bursts in this sample 
as well as a number of synthesised data sets we tested our code 
on. The 68% confidence interval of the Monte Carlo error anal- 
ysis should hence be treated as an approximate estimate of error. 



6. Conclusion 

We have developed a method whereby early time photometric 
nlR, optical and UV data of the afterglow can be used to estimate 
the photometric redshift of GRBs. This has advantages over the 
more usual methods of spectroscopic redshift determination, in 
that photometric observations do not require the source to be 
as bright, and over photometric estimates of the host galaxy as 
the afterglow is frequently much brighter than the host. In this 
implementation, we assume that GRB afterglow spectra are well 
represented by a power-law, and model the efi'ects of absorption 
due to the Lyman forest and host extinction. We use a genetic 
algorithm-based routine to simultaneously fit the parameters of 
interest, relatively unbiased by assumptions of the spectral index 



p. A. Cuiran et al.: Robust photometric redshift determinations of GRB afterglows at z ^ 2 





0.05 0,1 0.15 0.2 0.5 1 1.5 



Host extinction, Eg_y Spectral index, /? 

Fig. 5. Distribution plots for the Monte Carlo error analysis of GRB 050814. Redshift, host extinction and spectral index distributions are plotted 
for all trials of the Monte Carlo analysis, and overlaid on the re-binned results of the genetic algorithm optimum population from Figure 4 (in 
black). 



or host extinction. We introduce new nIR data for GRB 990510, 
which we have fitted along with the previously published data, 
to give new estimates of the host extinction and spectral index. 
We offer a new photometric redshift for GRB 050814, slightly 
higher than that previously suggested. 

The Monte Carlo error analysis, though computationally 
time consuming, is required as the method of error estimate 
based on the optimum population width of the genetic algo- 
rithm underestimates errors significantly. As the distribution of 
the best fit parameters obtained via Monte Carlo analysis are not 
Gaussian, caution is required when interpreting the nominal 68% 
errors of the average parameters. 
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